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ABSTRACT 


The  purpose  of  this  study  was  to  determine  analytically  the 
effects  of  water  vapor  on  the  ignition  delay  process  in  the 
burning  of  gaseous  hydrogen  in  a  supersonic  flowing  vitiated- 
airstream.  A  low  temperature  (950°K  «  1000°K)  regime  was  used 
so  the  reaction  rates  would  be  dominate  over  the  diffusion  rates. 

In  this  study  two  chemical  kinetic  computer  programs  were  used  to 
analytically  predict  the  ignition  delay  times  using  the  inputs 
from  seven  vitiate d-air  heater  experimental  investigative  con¬ 
ditions.  The  analytically  calculated  ignition  delay  times,  for 
the  seven  vitiated-air  and  seven  comparitive  clean-air  cases,  were 
then  compared  to  the  experimental  results  to  determine  the  effects 
of  <ater  vapor  on  the  ignition  delay  process. 

It  was  found  that  the  water  vapor  shortened  the  ignition  delay 
time  in  the  temperature  regime  of  this  investigation. 


A  STUDY  OF  NONEQUILIBRIUM  EFFECTS 
ON  IGNITION  DELAY  IN  A  SUPERSONIC 
AIRSTREAM 


Ie  INTRODUCTION 

Supersonic  ramjet  engines  (SCRAM JETS)  require  larga  mass 
flow  rates  of  air  to  produce  desirable  thrust.  Ground  test  facil¬ 
ities  which  could  provide  air  mass  flow  rates  at  temperatures, 
pressures,  and  velocities  necessary  for  meaningful  engine  develop¬ 
ment  are  virtually  non-existent. 

Due  to  floatation  and  erosion  problems  associated  with  high 
mass  flow  rates  in  clean  air  heaters,  vitiated  air  heaters  may  be 
the  only  solution  to  achieving  high  energy  fluids  necessary  for 
the  ground  development  of  SCRAM JET  engines. 

Combustion  type  vitiated  air  heaters  are  gas  generator  and 
nozzle  combinations  that  produce  large  mass  flow  rates  with  inherent 
high  temperatures  and  pressures  as  a  result  of  the  combustion  of 
propellants.  Careful  design  of  these  apparatuses  cm  duplicate 
flight  conditions  necessary  for  hypersonic  operation  of  a  SCRAMJET 
engine  in  a  ground  test  facility.  The  working  fluid  generated 
which  is  called  vitiated  air,  has  the  same  molecular  weight  and 
oxygen  content  as  clean  air.  Additional  species  such  as  water 
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vapor  and  carbon  dioxide  may  be  presen'  due  to  by-products  of  the 
combustion  process  and  are  considered  contaminants  when  compared 
with  real  air  composition*  If  the  effects  of  these  contaminates 
on  the  ignition  process  between  fuel  and  vitiated  air  could  be 
determined,  than  the  vitiated  air  heater  could  be  economically 
used  for  meaningful  engine  combustor  development  in  ground  test 
facilities. 

The  purpose  of  the  study  is  to  determine  analytically  the 
effects  of  water  vapor  contamination  on  the  ignition  process  of 
heated  gaseous  hydrogen  fuel  in  a  Mach  3  vitiated  airstrearn. 


Scope  of  Present  Investigation 


This  study  is  an  extension  of  the  work  performed  by  Major 
Renaldo  M.  Jensen  at  Purdue  University  as  partial  fulfillment  of 
the  requirements  for  the  Doctor  of  Philosophy  Degree  (Ref, 6  ). 
Jensen's  experimental  results  on  the  ignition  delay  phenomenon  are 
included  in  this  thesis  as  a  comparison  for  the  analytical  results 
determined  in  this  study. 

The  SIL  (Simplified  Ignition  Lag)  (Ref. 6  )  computer  program 
that  Jensen  used  to  analytically  predict  the  ignition  delay  lengths 
in  his  experimental  cases  did  not  allow  for  a  temperature  increase 
during  the  ignition  delay  period. 

It  was  felt  that  some  validity  in  the  comparison  of  his  ana¬ 
lytical  and  experimental  d  ta  was  impaired  by  neglecting  this  tem¬ 
perature  rise.  However,  in  this  study  the  ignition  delay  lengths 
of  Jensen’s  experimental  cases  will  be  recalculated  allowing  for 
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the  temperature  increase  during  the  ignition  delay  period  to  de¬ 
termine  whether  the  allowance  of  this  temperature  increase  im¬ 
proves  his  correlation  between  the  experimental  and  analytical 
results.  Additional  Ignition  delay  lengths  will  be  computed 
using  clean  air  rather  than  vitiated  air  for  the  same  experimental 
conditions,  to  determine  the  effects  of  the  water  vapor  on  the  ig¬ 
nition  delay  process. 

Review  of  Pertinent  Literature 

Jensen  (Ref. 6  ).  Jensen  designed  and  built  a  vitiated  air 
heater  (Fig.  1,  Pg. 4  >  that  reacted  a  nitrogen  (N2>  diluent,  nitro¬ 
gen  tetroxide  (HjOi,),  and  hydrazine  (N2Hj,)  to  generate  the  vitiated 
air.  The  use  of  nitrogen  (N2)  increased  the  molecular  weight  to 
that  cf  real  air.  The  water  vapor  inherently  produced  as  a  combus¬ 
tion  product  in  the  gas  generation  process  (19  percent  by  weight) 
replaced  an  equal  amount  of  nitrogen  in  the  vitiated  air. 

This  system  was  designed  to  duplicate  the  temperature  (1000°K), 
pressure  (14«7  PSIA),  and  velocity  conditions  (M  =  2,9)  of  air 
entering  the  combustor  of  a  SCRAMJET  engine  operating  at  a  flight 
mach  nuriber  of  6  at  an  altitude  of  80,000  teet.  This  simulation  of 
a  SCRAMJET  combustion  chamber  entrance  was  further  enhanced  by  using 
an  annular  injector  ax  the  exit  of  the  gas  generator  nozzle  for  the 
introduction  of  the  heated  hydrogen  gaseous  fuel  (590°  to  750°K) 
into  the  vitiated  air  stream.  Injection  was  parallel  to  the  flow 
in  order  to  minimize  the  recirculation  effects  between  the  concen¬ 
tric  streams  of  the  free  jet,  (Fig.  2,  Pg.6  ).  Once  ignition 
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occurred  In  this  vnconfined  flow  it  was  possible  to  physically 
Measure  tha  ignition  delay  phenomenon  captured  by  a  high  speed 
motion  picture  (Fastex)  camera  (5S0  frames  a  min.).  See  Fig  3, 
Pg.7  and  Fig.  4,  Pg.a  « 

The  S1L  computer  program  supplied  by  NASA  Lewis,  was  used  to 
analytically  predict  the  ignition  delay  lengths  for  his  experi¬ 
mental  cases.  It  was  found  that  the  ignition  delay  length  de¬ 
creased  with  increasing  hycrogen  inlet  temperature  and  that  the 
ignition  delay  lengths  were  comparatively  independev :  of  the  cqui- 
valence  ratio  (Figs.  5  and  6,  ?g.9  ).  The  results  also  showed 
that  the  presence  of  water  vapor  shortened  the  ignition  delay 
length  (Fig.  7,  Pg.ll). 

Snyder >  Robertson,  Zanders,  and  Skinner  (Ref. 8  ?. 

In  this  study  the  effects  of  water  vapor,  nitric  oxide,  ni¬ 
trogen  dioxide  and  ammonia  on  ignition  delay  were  studied  using 
shock  tube  techniques.  The  ranges  considered  were,  equivalence 
ratio  .5  to  1.0,  temperature  from  800°K  to  1100°K,  and  the  pres¬ 
sure  range  from  1.5  PSIA  to  130  PSIA. 

Their  results  showed  that  for  an  equivalence  ratio  of  1, 
pressure  of  30  PSIA,  and  water  vapor  content  of  20  percent,  that 
the  water  vapor  present  would  shorten  the  ignition  delay  if  the 
temperature  is  higher  than  970°K,  and  would  lengthen  the  ignition 
delay  if  the  teverature  is  lower  than  970°K  (Fig,  8,  Pg.13). 

Their  results  also  showed  that  nitric  oxide  (Fig,  8)  and  nitrogen 
dioxide  had  an  accelerating  effect  on  the  ignition  delay.  They 
also  indicated  that  ammonia  had  very  little  effect  on  the  ignition 
delay  in  this  low  temperature  range. 
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Figure  4.  IGNITION  DELAY  MEASUREMENT  USING  HIGH  SPEED 
FASTEN  CAMERA 
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Edelman  and  Spadaceini  (Ref. 4  )•  This  study  was  accom¬ 
plished  using  preheated  air.  The  air  was  preheated  by  burning 
oil  and  hydrogen  gas.  By  adding  make-up  oxygen  to  a  combustion 
heated  true  air,  (revitalized  air),  it  was  possible  to  obtain 
vitiated  air  with  the  correct  proportion  of  atmospheric  oxygen. 

They  concluded  from  their  experimental  and  analytical  results  that 
water  vapor  and  carbon  dioxide  could  either  increase  or  decrease 
the  ignition  delay  time  depending  on  the  initial  temperature,  pres¬ 
sure,  and  concentration  level. 

Approach 

The  ODFRC  (One  Dimensional  Finite  Rate  Chemistry)  computer 
program  will  be  used  to  analytically  determine  the  ignition  delay 
lengths  for  the  seven  vitiated  air  experimental  cases.  (Ref, 3  5  7). 
This  program  was  developed  oy  GASL  and  provided  for  this  research 
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program  by  NASA  Lewis  Research  Center.  It  differs  from  the  SIL 
Computer  program  used  by  Jensen  in  that  it  takes  into  account  the 
temperature  increase  during  the  induction  period. 

The  GCK  (General  Chemical  Kinetics)  computer  program  will  also 
be  used  to  determine  the  ignition  delay  lengths  for  the  seven  clean 
and  vitiated  air  cases.  This  program  was  developed  by  NASA  Lewis 
Research  Center  (Ref,  2  ).  The  GCK  program,  is  similar  but  slightly 
mere  sophisticated  than  the  ODFRC  computer  program  in  that  it  has 
the  added  ability  to  allow  for  up  to  thirty  reactions.  Also  unlike 
the  other  computer  programs  written  for  combustion  studies,  this 
program  internally  calculates  the  step  size  required  for  the  inte- 
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gration  of  the  rate  equations,  which  help  eliminate  the  time  prob¬ 
lem  that  exists  with  other  integration  schemes.  Also  this  program 
allows  for  no  reactions  which  arc  felt  to  be  impo  «:ant  in  ignition 
delay  calculations. 

To  determine  the  nonequilibrium  combustor  inlet  clean  and  viti¬ 
ated  air  species  concentration  used  as  the  data  input  to  the  ODFRC 
and  GCK  programs  the  FRRG  (Finite  Rate  Reacting  Gas)  computer  pro¬ 
gram  will  be  used.  This  computer  program  was  developed  by  TRW  Sys¬ 
tems  for  NASA  (MSC)  under  contract  NAS-4358  (Ref.  5).  It  was  also 
used  by  Jensen  to  calculate  the  inlet  species  concentrations  for 
use  in  the  SIL  computer  program. 

The  results  of  the  clean  and  vitiated  air  case  will  be  compared 
to  the  experimental  results  to  determine  the  effects  of  the  water 
vapor  on  the  ignition  delay  lengths. 


Assumptions 

In  the  S'  ■*.n  experimental  cases,  the  clean  and  vitiated  air 
entering  the  combustor  are  assumed  to  be  in  chemical  nonv-qui librium. 
It  is  also  assumed  that  the  airstream  entering  the  combustor  is  a 
homogeneous  mixture  since  the  computer  programs  are  not  designed  to 
handle  mixing  phenomena. 

The  assumptions  made  for  the  FRRG,  CDFRC  and  GCK  computer 
programs  are  as  follows: 

1.  The  flow  is  adiabatic. 

2.  The  species  behave  as  ideal  gases. 

3.  Viscous  effects  are  neglected. 
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II.  ANALYTICAL  PROGRAM 


The  author  used  the  FRRG  program  to  analytically  determine 
the  nonequilibrium  composition,  temperature,  pressure,  and 
velocity  characteristics  of  the  working  fluid  expanding  in  the 
gas  generator.  The  computed  outputs  were  then  used  as  the  input 
data  to  the  ODFRC  and  GCK  programs  to  analytically  predict  the 
ignition  delays  in  the  free  jet  which  simulates  the  ignition  pro¬ 
cess  in  a  SCRAM JET  combustor. 

Description  of  Computer  Programs 

The  computer  programs  indicated  in  the  above  section  are 
discussed  in  the  following  paragraphs  under  their  individual 
headings.  Included  in  these  discussions  are: 

1.  Program  description 

2.  Input  requirements 

3.  Output 

4.  Limitations 

Finite  Rate  Reacting  Gas  (FRRG)  Program.  The  FRRG  program 
was  developed  ’jy  TRW  Systems  for  NASA  (MSC)  under  contract  NAS-435 C. 
This  program  solves  the  equations  for  invicid  one  dimensional  flow 
of  a  reacting  gas  mixture  (Ref.  5).  To  minimize  computer  time  during 
the  integrations  of  the  chemical  relaxation  equations,  the  program 
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uses  an  implicit  integration  scheme  for  near  equilibrium  flow, 
and  explicit  integration  scheme  for  near  equilibrium  flow,  and 
explicit  integration  for  flow  distant  from  equilibrium. 

Required  data  inputs  are  the  nozzle  geometry,  propellant 
data,  and  reaction  rate  information.  The  reaction  rate  expres¬ 
sions  and  rate  constants  used  in  this  program  are  listed  in  Table 
1,  pg.  18  .  Only  the  forward  rates  are  required  input  since  the 
reverse  reaction  rates  are  internally  calculated  from  the  forward 
reaction  rate  and  the  equilibrium  constant  by  the  following  ex¬ 


pression! 


•  =  kf<J)/K  (j) 


Also  when  the  reaction  equations  used  were  the  same  as  in 
the  FRRG,  ODFRC,  and  GCK  programs,  identical  input  reaction  rate 
data  was  used  to  maintain  continuity  in  the  results. 

This  program  has  the  option  of  performing  either  equilibrium 
or  finite  rate  calculations.  Output  from  this  program  are  tem¬ 
perature,  pressure,  maoh  number,  velocity,  density,  enthalpy,  gas 
molecular  weight,  heat  capacity,  ratio  of  specific  heats,  and 
species  concentrations.  The  finite  rate  calculation  was  employed 
since  this  option  more  closely  duplicated  the  experimental  con¬ 
ditions. 

Output  species  concentrations  are  printed  in  mole  fractions 
and  mass  fractions.  Mass  fraction  units  were  chosen  since  the 
ODFRC  program  requires  that  the  species  concentrations  be  input 
in  terms  of  mass  fractions. 

One  Dimensional  Finite  Rate  Chemistry  (ODFRC)  Program. 

The  ODFRC  program,  supplies  by  RASA  Lewis  Research  Center,  solves 
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for  the  rate  of  growth  of  the  free  radicals  as  a  function  of  time 
by  using  a  standard  matrix  technique.  These  rates  are  then  inte¬ 
grated  to  give  the  species  concentrations  at  that  instant  of  time. 
Unlike  the  SIL  program  used  by  Jensen  (Ref.  6)(  this  program  al¬ 
lows  for  the  temperature  increase  during  the  ignition  delay  period 
caused  by  certain  reactions  being  slightly  more  exothermic  than 
the  other  reactions  are  endothermic* 

Curing  the  development  of  the  ODFRC  program,  modifications 
were  made  by  NASA  Lewis  Research  Center  to  the  original  program 
to  allow  the  computation  to  run  though  the  steep  temperature  rise 
following  ignition. 

The  program  has  the  option  of  using  pressure  vs,  distance, 
area  vs.  distance,  or  radius  vs.  distance  for  combustor  geometry 
conditions.  The  pressure  vs.  distance  option  was  chosen  to  dupli¬ 
cate  the  constant  pressure  experimental  conditions  of  a  free  jet 
burning  in  the  atmosphere.  Other  inputs  include  initial  temper¬ 
ature,  pressure,  velocity,  area,  distance,  time,  species  concen¬ 
trations,  and  reaction  rate  constants.  Reactions,  rates,  and 
constants  used  are  listed  in  table  2,  page  21  .  Species  considered 
were  H,  0,  HjO,  H2,  O2,  OH,  N2,  HO2,  and  H2O2, 

Output  from  this  program  includes  time,  temperature,  density, 
pressure,  enthalpy,  velocity,  area,  distance,  and  species  concen¬ 
trations  for  the  nine  species  listed  above.  The  author  modified 
the  program  to  include  a  calcomp  plotting  subroutine  to  plot  0  H 
concentration  vs,  time,  H  concentration  vs,  time,  and  temperature 
vs,  time.  The  plots  which  are  used  to  determine  the  ignition  de¬ 
lay,  will  be  discussed  in  a  later  section. 
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General  Chemical  Kinetics  (GCK)  Program.  The  GCK  computer 
program  was  supplied  by  NASA  Lewis  Research  Center  for  this  study. 
Unlike  the  ODFRC  Program,  the  GCK  program  allows  for  reactions 
involving  NO  which  possibly  cculd  have  an  accelerating  effect  on 
the  ignition  delay. 

The  program  solves  the  reaction  equations  by  :ans  of  an  im¬ 
plicit  integration  technique.  This  implicit  solution  scheme  was 
combined  with  a  new  step-size  optimization  method.  This  new  inte¬ 
gration  technique  which  is  efficient  for  both  ve*y  slow  and  very 
fast  (near-equilibrium)  chemical  reactions  (Ref. 2  )•  Step-size 
•control  decreases  the  computation  time  required. 

Input  conditions  for  this  program  allow  for  30  reactions. 

The  reactions  and  reaction  rate  constants  used  in  this  program 
are  listed  in  table  3,  pg.23  .  Reverse  reaction  rates  are  calcu¬ 
lated  internally  from  the  relation: 


Kj,(J)  =  kf(J>  /KC(J) 


In  addition  to  the  chemical  species  which  are  required  by  the 
given  reactions,  the  user  can  include  one  or  more  inert  species, 
but  the  total  number  of  species  considered  cannot  exceed  25  (Ref. 2). 
For  each  dissociation-recombination  reaction,  the  user  can  specify 
a  set  of  third-body  efficiencies  (Ref. 2  ).  To  maintain  continuity 
between  the  computer  programs,  the  same  efficiencies  were  used  in 
the  GCK  program  as  in  the  ODFRC  program  where  possible.  If  the 
third-body  efficiencies  are  not  specified  they  are  assumed  to  have 
a  value  of  one. 
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Initial  conditions  required  are: 

1«  Time  and  axial  position. 

2.  Velocity  (or  mach  number). 

3.  Density  (or  pressure). 

4.  Temperature. 

5.  Species  c  mcentrations 

Also  pressure  or  area  profile  must  be  input  if  the  problem  under 
consideration  is  a  flow  problem.  (Ref. 2) 

Outputs  from  this  program  are  temperature,  time  pressure, 
velocity,  mach  number  and  species  concentrations.  The  output  can 
be  in  CGS,  TPS,  or  SI  units*  The  CGS  option  was  chosen  to  maintain 
continuity  with  the  results  of  the  other  programs. 
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Analytical  Determination  of  the  Ignition  Delay  Time. 


The  PRRG  program  was  vised  to  calculate  the  species  concen¬ 
trations,  velocity,  temperature  and  pressure  conditions  at  the 
exit  of  the  vitiated  air  heater.  The  exit  of  the  vitiated  air 
heater  is  considered  the  inlet  of  the  combustor  in  this  study. 
Tha  species  mass  fractions  were  recalculated  allowing  for  the  H2 
gas  that  is  supplied  at  the  combustor  inlet.  These  species  mass 
fractions,  mixing  cup  temperature,  and  velocity  conditions  are 
used  as  the  input  to  the  ODFRC  and  GCK  programs.  The  C0FRC  and 
GCK  calculate  the  temperature,  velocity,  and  species  concentra¬ 
tion  as  a  function  of  time.  See  Tig,  9,  pg.  24. 

The  end  of  the  ignition  delay  period  is  characterized  by  a 
rapid  increase  in  termperature  (Fig.  10,  Pg,26).  During  the 
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TABLE  3 


REACTION  RATES  FOR  GCK  PROGRAM 


REACTION 


H  ♦  02  =  OH  =  0 
0  ♦  H2  *  OH  ♦  H 
OH  =  H20+  H 


0  +  HjO  =  OH  +  OH 


h2+  ho2  =  h2o2+  H 

H2  ♦  02  =  OH  ♦  OH 

H  +  HO„  =  OH  +  Oh 
2 

8  2H0  s  HO  *  0^ 

2  2  2  2 

H  +  H202=  H20  +  OH 

10  OH  +  H202  =  H20  +  H02 

11  M  +  H  +  OH  =  HjO  +  H 

12  M  +  K  ♦  02  =  H02  +  H 

13  MtH+H  s  Hj  t  K 

14  Kv02+Os0  +  M 

15  M  +  H202  =  OH  ♦  OH 

16  M  +  N  ♦  N  s  N2  +  M 

17  M  +  N  +  0  =  NO  +  M 


A  (J) 
(gtn.cm) 


1.25  X  1014 

1.7  X  1013 
2.19  X  1013 

5.75  X  1013 
9.6  X  1012 
1.0  X  1013 
7,0  X  1013 

1.8  X  1012 
3.1  X  1014 
1.0  X  1013 
7.5  X  10ie 
1.59  X  1015 
1.0  X  1018 

2.75  X  1019 
1.17  X  lOl? 

2.8  X  i0l7 
5,4  X  1016 


B(J) 

(kcal/mole)  N(J) 


s. 
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Figure  9.  PROGRAM  APPLICATION 


GAM/ME/72-1 


Ignition  delay  region  the  net  results  of  the  reactions  are  slight¬ 
ly  exothermic*  This  accounts  for  the  slight  temperature  increase 
during  the  ignition  delay  period.  At  the  point  where  the  temper¬ 
ature  increases  rapidly,  the  recombination  reactions  are  the  con¬ 
trolling  reactions.  Since  the  recombination  reactions  are  exother¬ 
mic,  the  temperature  increases  rapidly  (Fig,  10,  26). 

The  following  criteria  were  used  to  determine  the  time  at 
which  the  ignition  delay  period  ended. 

1.  Intersection  of  the  linear  slope  temperature  line 
in  the  ignition  delay  zone,  with  the  linear  slope  tem¬ 
perature  line  of  the  reaction  zone.  (Fig,  11,  Pg.  27) 

2.  Deviation  from  the  linear  portion  of  the  OH  vs. 
time  curve  (Fig,  12,  Pg.  28) 

3.  Deviation  from  the  linear  portion  of  the  H  vs. 
time  curve  (Fig,  13,  Pg.  29) 

4.  One  half  of  the  initial  H2  concentration. 

All  four  methods  yield  essentially  identical  delay  times  so  all 
methods  appear  applicable. 
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Figure  10.  TEMPERATURE  HISTORY  DURING  IGNITION 


Figure  12.  GRAPHICAL  REPRESENTA f I ON  OF  IGNITION  DEI, AY  USING  OH  CRITERIA 
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III,  DATA  REDUCTION 


Analytical  results  were  obtained  by  using  the  OH  concen¬ 
tration,  H  concentration,  concentration,  and  the  temperature 
slope  criteria  as  discussed  in  section  II.  Due  to  the  data  simi¬ 
larity,  only  the  OH  concentration  criterion  results  will  be  dis¬ 
cussed  at  length.  Graphical  data  based  on  the  remaining  criteria 
are  given  in  the  appendix.  Discussion  of  OH  data  results  is  con¬ 
sidered  under  the  following  headings: 

1.  ODFRC  program  results 

2.  GCK  program  results 

ODFRC  Program  Results 

The  ODFRC  results  obtained  by  holding  the  equivalence  ratio 
constant  (1/  and  varying  the  inlet  hydrogen  gas  temperature  are 
shown  in  Fig.  14,  Pg.  31.  The  results  obtained  by  holding  the 
inlet  hydrogen  tenperature  constant  (710°K)  and  varying  the  equi¬ 
valence  ratio  are  shown  in  Fig.  15,  Pg.  32.  Note  that  the  data 
points  for  the  analytically  calculated  clean  and  vitiated-air 
cases  are  from  50%  to  100%  lower  than  the  experimental  values. 
Results  for  the  seven  cases  were  analytically  recalculated  using 
a  lower  temperature  to  see  if  the  ignition  delay  times  could  be 
increased.  All  the  analytical  data  failed  to  indicate  ignition 
at  thsse  1c  ix  temperatures. 


400  450  500  550  600  650  700  750  800 

HYDROGEN  TEMPERATURE  (°K) 

Figure  14.  ODFRC  IGIIITIOH  DELAY  HISTORY  USING  OH  CRITERION 
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.  Results  from  the  ODFRC  program,  when  comparing  clean  and 
vitiated-air  results,  show  that  the  presence  of  water  vapor  short¬ 
ened  tne  ignition  delay  time  (Figs.  14  and  15).  However,  due  to 
the  deviation  between  the  analytical  and  experimental  data  points 
and  the  added  sophistication,  the  GCK  computer  program  was  used 
in  an  attempt  to  analytically  predict  a  mori-  favorable  comparison 
between  the  analytically  calculated  vitiated-air  and  the  seven 
experimental  cases. 


GCK  Program  Results 


Fig.  16,  Pg.  35  and  Fig.  17,  Pg.  36  show  the  computed  results 
obtained  from  the  GCK  program.  These  values  compare  favorably 
with  the  results  of  the  ODFRC  program,  but  were  still  much  lower 
than  the  experimental  points.  The  mixing  cun  termperature,  which 
simulates  the  rate  of  diffusion  of  the  !I  into  the  working  fluid, 
was  decreased  to  the  value  that  would  compare  with  the  40%  mixing 
conditions  assumed  by  Jensen  (Ref. 6  )•  The  results  obtained  (Fig. 

18  and  19)  were  closer  to  the  experimental  results  than  the  values 
obtained  from  the  ODFRC  program,  but  the  slope  of  the  vitiatcd-air 
ignition  delay  curve  was  much  less  than  the  slope  of  the  experi¬ 
mental  curve.  A  possible  explanation  is  that  diffusion  rates  are 
much  slower  at  low  temperatures  than  at  high  temperatures  (Ref.  1). 
Decreasing  the  diffusion  rate  at  the  lower  temperatures  would  in¬ 
crease  the  slope  of  the  vitiated-air  ignition  delay  values  on  the 
constant  equivalence  ratio  curve  (Fig,  18)  and  increase  the  ignition 
delay  values  for  the  vitiatcd-air  cases  on  the  constant  temperature 
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*•  (Fig.  *>.  «-  vltlated-air  values  —  ■«-  *  —  *° 

the  experimental  «M  These  results  are  essentially  the  sane 
.,  the  analytic  results  Jensen  (Ref.  6)  ehtained  using  the  «L 
conputer  pregran  with  the  h0\  nixing  assumption. 


EXPERIMENTAL 
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Figure  17.  GCK  IGNITION  DELAY  HISTORY  USING  OH  CRITERION 


Figure  19.  GCK  IGNITION  DELAY  HISTORY  USING  OH  CRITERION  (40%  MIXING) 


GAM/ME/72-1 


IV.  CONCLUSIONS 

Analytical  calculations  ware  made  for  seven  vitiated  an.* 
clean-air  investigative  cases.  These  results  were  compared  to 
the  seven  low  temperature  supersonic  experimental  data  points 
obtained  by  Jensen.  Based  on  the  results  in  section  III,  the 
following  conclusion  was  made: 

Ignition  delay  time  is  decreased  if  water  vapor 
is  present  at  the  inlet  to  the  combustor.  Water  vapor 
acts  as  a  third  body  which  promotes  recombination  of  the 
free  radical  species  causing  the  exothermic  termination 
reactions  to  become  dominent. 
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V.  RECOMMENDATIONS 

A*  a  result  of  this  study,  the  following  recommendations 
are  offered: 

1*  Obtain  experimental  inputs  and  results  from  vitiated- 
air  heater  ignition  delay  studies  for  a  wider  temperature 
range  and  analytically  predict  the  ignition  delay  times. 
Results  of  this  study  would  determine  the  temperature 
regime  in  which  the  water  vapor  shor*-'"—  *  Onitio«  delay 
times. 

2,  Determine  the  rate  at  which  the  • '  into  the 

wording  fluid  for  different  temperature  0es.  Results 
from  this  study  could  be  used  to  more  accurately  predict 
ignition  delay  times. 

3.  Develop  a  kinetic  scheme  which  would  allow  for  more 
accurate  evaluations  of  the  effects  of  nitric  oxide  (NO), 
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Figure  21.  ODFRC  IGNITION  DELAY  HISTORY  ^SIHG  1/2  H-  CRITERION 
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Figure  24.  ODFRC  IGNITION  DE1AV  HISTORY  USING  1/2  H  CRITERION 
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Figure  25.  ODFRC  IGNITION  DELAY  HISTORY  USING  TEMPERATURE  SLOPE  CRITERION 
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